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1 Introduction 

Nowadays the considerable advance [IHS] in the rigorous derivation of the Boltzmann kinetic 
equation in the Boltzmann- Grad scaling limit |6l|7] is well known. The lack of similar progress 
for the Enskog kinetic equation [8HT2] suggested by D. Enskog [13] as a genelalization of the 
Boltzmann equation on dense gases, is conditioned by a priori stated collision integral of this 
kinetic equation. Thus, the rigorous derivation of the Enskog-type equations from dynamics of 
hard spheres [1], [H] remains an open problem so far. 

In this paper we consider the problem of potentialities inherent in the description of the 
evolution of states of a hard sphere system in terms of a one-particle distribution function. In 
view of the generic initial data of hard spheres are determined by initial correlations connected 
with the forbidden configurations of particles, one of the problems consists in the derivation of 
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the kinetic equations from underlaying dynamics in case of presence of correlations at initial 
time. We establish that in fact, if initial data is completely defined by a one-particle distribution 
function, then all possible states of infinitely many hard spheres at arbitrary moment of time 
can be described within the framework of a one-particle distribution function by the kinetic 
equation without any approximations. 

We briefly outline the main results and structure of the paper. In section 2 we formulate 
some preliminary facts about dynamics of finitely many hard spheres and state the main result. 
In sections 3 and 4 the main results related to the origin of the Enskog kinetic evolution are 
proved. On the basis of the kinetic cluster expansions of cumulants of groups of operators 
of finitely many hard spheres which give the possibility to construct the kinetic equations in 
case of presence of correlations at initial time, in section 3 we derive the generalized Enskog 
equation. Then in section 4 we establish that all possible states are described by a sequence 
of explicitly defined functionals of a solution of constructed kinetic equation. Moreover, in 
section 5 we state the Markovian generalized Enskog equation and establish its links with the 
Enskog- type kinetic equations [8], [TSHIZ]. In section 6 the existence theorem is proved for 
the abstract initial- value problem of the generalized Enskog equation in the space of integrable 
functions. Finally, in section 7 we conclude with some observations and perspectives for future 
research. 

2 On the origin of the Enskog kinetic evolution 

2.1 Preliminaries: dynamics of finitely many hard spheres 

We consider a system of a non-fixed, i.e. arbitrary but finite, number of particles of a unit mass 
(the formalism of nonequilibrium grand canonical ensemble |I]), interacting as hard spheres 
with a diameter a > 0. Every particle is characterized by the phase coordinates {qi,Pi) = 
Xi G M^ X M^, i > 1. For configurations of such a system the following inequalities hold: 
\qi — Q'jl > cr, « 7^ j > 1, i.e. the set W„ = {(q'l, . . . , gn) ^ IR^^I |(?i ~ (ij\ < cr for at least one pair 
{i,j) '■ i 7^ j G (1, . . . , r?,)} is the set of forbidden configurations of n hard spheres. 

The evolution of all possible states is described by the sequence of marginal (s-particle) dis- 
tribution functions Fs{t, Xi, . . . , Xg), s > 1, that satisfy the initial-value problem of the BBGKY 
hierarchy [1] 

r\ S p 

-F,{t) = C*Mt) + J2 / dx,^i^l,{t, s + l)F,+i(t), (1) 

F,{t)\t=o = F^, s>l. (2) 

If t > 0, in hierarchy of evolution equations ([1]) the operator C* is defined by the Poisson 
bracket of noninteracting particles with the corresponding boundary conditions on dWg |1]: 

S p. 

^*sFs{t) = -J2(Pi^ jr-)\aw,Fs{t,xi, . . . , Xs), (3) 



Enskog Equation 



where (77, {pi — ps+i)) = J2a=iV°'iPf ~ P^+i) is a scalar product. Operator ([3]) is the generator 
of the Liouville equation for states and it is an adjoint operator to the generator Cs of the 
Liouville equation for observables [Ij. The operator C*^^{i, s + 1) is defined by the expression 



J2 [ dxs+iCUi, s + l)F,+i(t) =a^J2 [ 



dps+idr] {r], {pi - Ps+i)) x (4) 



X 



[Fs+iit, Xi,...,qi,p*,...,Xs,qi- (tv,P*s+i) - Fs+i{t, Xu...,Xs,qi + ar],ps+i)) , 



where S^ = {r/ G M^| |?7| = 1, (77, {pi — Ps+i)) > 0} and the momenta p*, p^^^ are defined by the 
following equalities 

Pi =Pi-V {V, {Pi - Ps+i)) , (5) 

pI+1 = Ps+1 + V {V, {Pi - Ps+i)) ■ 

If t < 0, the generator of the BBGKY hierarchy is defined by the corresponding expression [1]. 
We will consider initial data (E]) satisfying the chaos condition [1], i.e. statistically indepen- 
dent hard spheres 

s 

Fs{t)\t=o = X{Fl{x^)X^^s\Y^,, (6) 

i=l 

where Afjjss^iyy^ = Xsi.Qi', . • . , ^s) is a characteristic function of allowed configurations M^* \ W^ of 
s hard spheres. Initial data ([6]) is intrinsic for the kinetic description of many-particle systems 
because in this case all possible states are characterized by means of a one-particle marginal 
distribution function. 

To construct a solution of initial- value problem (II])- (E]) we adduce some preliminaries about 
hard sphere dynamics. Let L\ = L^(]R^" x (M^" \ W„)) be the space of integrable functions 
that are symmetric with respect to the permutations of the arguments Xi, . . . ,a:„, equal to 
zero on the set of forbidden configurations W„ and equipped with the norm: ||/n||Li(R3nxK3n) = 
/ dxi . . . dxn\fn{xi, . . . , Xn)\- We denote by L;'^ q C L^ the subspace of continuously differen- 
tiable functions with compact supports. 

On a set of measurable functions the following one-parameter mapping: M 9 t t— )■ Sn{—t)fn 
is defined by 

['^n[~t) JnjyXl, . . . , Xn) = ^n[~t, 1, . . . , n)Jn[Xi, . . . , X^) = (7j 

'/„(Xi(-t,Xi,...,X„),...,X„(-t,Xi,...,X„)), if(xi,...,x„)G(M3"x(M3n\wj), 
0, ii{qi,...,qn)eWn, 

where Xj(— t) is a phase trajectory of ith particle constructed in [1]. We note that the phase 
trajectories of a hard sphere system are determined almost everywhere on the phase space 
R^"- X (M^" \ Wn) beyond the set M° of the zero Lebesgue measure [Hill]. The set M° consists 
of the phase space points which are specified such initial data that during the evolution generate 
multiple collisions, i.e. collisions of more than two particles, more than one two-particle collision 
at the same instant and infinite number of collisions within a finite time interval. 
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On the space L^ mapping l^^ defines the isometric strong continuous group of operators P, 
and on L^ Q C L^ the infinitesimal generator of group (ITj) is operator (^ which is the generator 
of the Liouville equation for hard spheres [1]. 

If F° G L^, then a nonperturbative solution of initial- value problem (II]), ([6]) is a sequence of 
distribution functions Fs(t,Xi, . . . , Xg), s > 1, represented by the following series [1], [H] 



Fs{t,xi,...,Xs) = (8) 

CO _ „ s+n 

= J2- I dxs+i . . . dxs+n^i+n{-t, {Y}, X\Y)Y\ F°(x,)A']R3(.+n)\ 



„^0^! J(R3xR3)n .^, 



where the evolution operator 2ti+„(— t) is the {n + l)t/i-order cumulant of groups of operators 
d?]) defined by the expansion 

2ti+„(-t,{y},x\r)= J2 (-i)i''-^(|p|-i)! n %x.)iH,^™), (9) 

P:{{Y},X\Y)=[j,X, X.CP 

and the following notations are used: {Y} is a set consisting of one element Y = (1, . . . , s), 
i.e. \{Y}\ = 1, ^p is a sum over all possible partitions P of the set {{Y},X \Y) = {{Y}, s + 
1, . . . , s + n) into |P| nonempty mutually disjoint subsets Xi G {{Y}, X \Y), the mapping 9 
is the declasterization mapping defined by the formula: 9{{Y},X \Y) = X. The simplest 
cumulants ([9]) have the form 

^,{-t,{Y}) = Ss{-t,Y), 

2l2(-t, {Y},s + l) = Ss+i{-t, y, 5 + 1) - S,i-t, Y)Sii-t, s + l). 

If ||-F{'||j;^i(K3xR3) < e^^, series (^ converges in the norm of the space L\ for arbitrary t G M. 
We remark that this condition means the condition on average number of particles. If we 
renormalize the distribution functions Fs{t) = -^Fg^t), where - is the density of particles (the 
average number of particles in a unit volume), we obtain solution expansion (|8]) over powers of 
the density ^ which converges provided that: ^ < e~"'^||-F°||~wjg3^jg3.. 

We note that nonperturbative solution ([8]) of the BBGKY hierarchy ([T]) is transformed to 
the form of the perturbation (iteration) series as a result of applying of analogs of the Duhamel 
equation [27] to cumulants of groups of operators Qj. 

In case of initial data F° G L^ g sequence of functions ([8]) is a strong solution of initial- value 
problem ([I]), ([6]) and for arbitrary initial data F^ G L\ it is a weak solution [18j. 

In the next section we describe the evolution of states within the framework of the kinetic 
theory, i.e. in terms of a one-particle marginal distribution function. 

2.2 The kinetic evolution of hard spheres: main results 

In view of the fact that initial data iQ is completely determined by a one-particle marginal 
distribution function on allowed configurations, the states given in terms of the sequence 
F{t) = {Fi{t, xi), . . . ,Fs{t,xi, . . . , Xs), . . .) of marginal distribution functions (IE]) can be de- 
scribed within the framework of the sequence F{t \ Fi{t)) = (Fi(t, xi), F2{t, xi, X2 | Fi(t)), . . . , 
Fs{t,xi, . . . ,Xs I Fi(t)), . . .) of the marginal functional of the state Fs{t,xi, . . . ,Xs \ Fi{t)), s > 
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2 which are exphcitly defined with respect to the solution Fi{t,xi) of the kinetic equation. 
We refer to such kinetic equation of a hard sphere system as the generahzed Enskog kinetic 
equation. 

Indeed for t > the one-particle distribution function ([8]) is a solution of the following 
initial-value problem 

-Fi(t,xi) = -(pi, — )Fi(t,xi) + (10) 

+a^ V] — / dp2dr] / dx^ . . . dxn+2{v, {Pi - P2)) x 

n=0 J^^^x^l J(R3xR3)n 

n+2 

x{^^+^{t,{l*,2*_},3,...,n + 2)F,{t,q^,pl)Fi{t,q,-ar],p;)l[F^{t,x,)- 

1=3 

n+2 

-5Ji+„(t, {1, 2+}, 3, . . . , n + 2)Fi(t, xi)Fi(f, gi + ar7,p2) n ^i(^' ^i)) ' 

j=3 

F,{t,x,)\t=o = F^ixi). (11) 

In this evolution equation we use notations from definition (jTj) adopted to the conventional 
notation of the Enskog collision integral: indices (l'',2!|_) denote that the evolution operator 
2Ji+„(t) acts on the corresponding phase points (^1,^1) and {qi ± ar],pl), and in general case 
the (n -|- l)t/i-order evolution operator QJi+„(t), n > 0, is defined by the expansion {\Y\ = 2 in 
case of equation (TTOj) ) 

n n n-mi-...-m^,_i 

ffl,«(Mn,x\r).„!$:(-i)'E- E („_„,_,.._,„,), >' (12) 

fc=0 -mi=l mj.=l 

x2li+n_mi-...-mfc(t, {^},s + 1,. . .,s + n-mi - ... - mfc) x 

/fc m,j ""n — mx-...-mj+s — 1 s-\-n — mi~. 

HE E n 

J=i ki=o fci_„^_..._„^+,=o v=i 
x2ii+fc^^_ , -ki ^ (t,ij,s + n-mi-...-mj + l + 



/fc rrij "'n—mi-...-mj+a—i s+n—m\~...~mj 

^n^- ^ TT 1^^ -X 

7=1 uJ—n U _n «9=1 V'^n-mi-...-mw+s+1-i,- "'n-7ni-...-mw+s+2-i,- /• 



n — 77T.1 — ... — m„'+s + l — i^- n — 7?^-l — ...— 7n,,+s + 2 — i„- 



where we mean fcj = rrij, fc^_^j_ ..^.^s+i = 0, and by the operator 2li+„(t) we denote the 
(1 + n)t/i-order scattering cumulant 

s+n 

2ti+„(t, {y}, X\Y) = 2li+„(-t, {F}, X \ y) J,+„(X) J] 2li(t, 2). (13) 

j=l 

In definition ( 1T3|) the operator 2li+„(— t) is the (1 + n)t/i-order cumulant (|9]) of groups ([7]) and 
the operator Jg+n is defined by the formula 

^s+n{X)fs+n = '^R3(s+n)\Ws+„/s+ni (14) 
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where '^Rscs+nj^w^^^ is a characteristic function of allowed configurations of s + n hard spheres. 
In section 4 we will give an illustration of evolution operators ( TT2l) . 
If t < 0, the generalized Enskog equation has the corresponding form 

— Fi(t,xi) = -(pi,— )Fi(t,xi)+ (15) 

+a^ V" — / dp2dr] / dx^ . . . dxn+2{'n, (Pi - P2)) x 

n=0 "']R^X§2_ J(K3x]R3)n 

n+2 

X (2Ji+„(t, {1, 2_}, 3, . . . , n + 2)Fi(t, qi,pi)F,{t, gi - (T?7,p2) JJ ^i(^' ^0 - 

i=3 

n+2 

-5Ji+„(t,{r,2;},3,...,n + 2)Fi(t,gi,pt)i^i(t,gi + a77,p;)J]Fi(t,x,)). 



i=3 



The marginal functional of the state Fs(t,xi, . . . ,Xs \ Fi{t)) , s > 2, are represented by the 
following expansions over products of the one-particle distribution function Fi (t) 

Fs{t,xi,...,Xs\Fi{t))= (16) 

00 ^ „ s+n 

= V^ / dXs+l...dXs+n'^l+n{t,{Y},X\Y)T\Fi{t,Xi), 

n=0 ''■ Ji^'x^')" ^=l 

where the {n + l)t/i-order evolution operator 5Ji+„(t), n > 0, is defined by expansion flT2|) . 
If ||-Fi(t)||ii(u3xiR3) < e~^^^~^'^\ series (IT6|) converges in the norm of the space L] for arbitrary 
teR. 

We remark that in terms of marginal functional of the state flT6|) the collision integral of the 
generalized Enskog equation is represented by the expansion 

Igee = ctM dp2dr] [r], (pi - P2)) x 

x(F2(f,gi,Pi,gi -(Tri,p* \ Fi(t)) -F2(t,Xi,gi + ar],p2 \ Fiit))). 

In section 5 we establish the relationship of the generalized Enskog equation (I10j) . (ll5p with 
the Markovian generalized Enskog equation and the revised Enskog equation [TSIIIS]- Here 
we only point out that first term of the collision integral expansion in kinetic equation fITOl) is 
exactly the collision integral of the Boltzmann- Enskog equation. 

Thus, under the condition: ||-Fi°||li(r3xr3) < e~'^^'*"'"^^(l + e~('^''+'^))~^, initial-value problem 
of the BBGKY hierarchy ([I]),([n]) is equivalent to initial- value problem of the generalized En- 
skog equation flTIl - flTTl) and the sequence of marginal functional of the state Fs{t,xi, . . . ,Xs \ 
Fi{t)), s > 2, defined by expansions (TT6l) . The proof of this statement is the subject of the next 
two sections. The existence theorem of the Cauchy problem (fTO!) - (fTT!) in the space of integrable 
functions is stated in section 6. 

Finally, it should be noted that the possibility to describe the evolution of all possible states 
of hard spheres in an equivalent way by the restated Cauchy problem for the generalized En- 
skog equation and by a sequence of the marginal functionals of the state as compared with the 



Enskog Equation 8 

corresponding Cauchy problem of the BBGKY hierarchy is an inherent property of the descrip- 
tion of many-particle systems within the framework of the formalism of nonequilibrium grand 
canonical ensemble which adopted to the description of infinite-particle systems in suitable 
functional spaces [1]. 

2.3 A one-dimensional hard sphere system 

We consider a one- dimensional system of hard spheres (hard rods of the length cr > 0). In view 
of that in this case for t > operator (jl]) is defined by the expression 



V / dxs+iCl^S, s + l)F,+i(t) = 

V" / dPP{F,+i{t,Xi, . . . ,qi,pi - P, . . . ,Xs,qi- cr,p^) - Fs+i{t,Xi, . . . ,Xs,q^ - a,p^ + P) + 

+F,+i(t, Xi,...,qi,Pi + P,...,Xs,qi + cr,Pi) - Fs+i{t, Xi, . . . ,Xs,qi + cr,Pi - P)), 

then for t > the generalized Enskog kinetic equation for a one-dimensional hard sphere system 
has the form 

— Fi(t,xi) = -pi— Fi(t,xi) + 

+ y2^ dPP{F2{t, quPi -P,qi- (y.Pi I F,{t)) - F2it, q^,p,, q,-a,p, + P\ F^{t)) + 

Z^r.^- Jo 



+F2it,quPi + P,qi + a,pi \ Fi{t)) ~ F2it,qi,pi,qi + a,pi - P \ Fi(t))), 

where the marginal functional of the state F2(t,xi,X2 \ Fi{t)) is defined by formula flT6l) . The 
series in right-hand side of this equation converges if ||-Fi(t)||Li(i8x]R) < e~^ (see section 4). For 
t < 0, the structure of the collision integral of the generalized Enskog equation induces by the 
structure of the corresponding term of the BBGKY hierarchy for a one-dimensional hard sphere 
system 



3 The generalized Enskog kinetic equation 

3.1 AuxiUary facts: kinetic cluster expansions 

To reformulate the Cauchy problem ([I]), ([6]) as the new Cauchy problem for a one-particle 
distribution function governed by the kinetic equation and the sequence of explicitly defined 
marginal functionals of the state which determined by a solution of such Cauchy problem, we 
introduce the following cluster expansions of cumulants (^ of groups of operators ([7]) 

2li+„(-t, {Y}, X \ Y)3s+„iX) = (17) 

n , 

V — \ 72' 

= 2^ (n-k)\k ! ^i+"-^i(^' ^^^' s + 1, • • • , s + n - A;i) X 
fci=0 ^ ^'' ^' 



Enskog Equation 



fcl , I fen-fc]^+s-l . 



^ k2\{ki - k2)\ " ' ^-^ kn-k^+sKkn-ki+s-l - kn-ki- 
/C2=U k„-ki+s='J 



11 '^i+k„.k^+s+i-^-k„-k^+,+2-i{-'t, ^,s + n-kl + l + ks+n-ki+2-i, ■ ■ 



s+n—ki 
X 

...,s + n-ki + ks+n-k^+i-ipi+k„_^^+,+,_i-k„_^^+,+2.M, s + n-ki + l + 

+ "^s+n-fci+2-i) ■ ■ ■ , S -\- n — ki + ^s+n-fei+l-ij) "^ ^ 0, 



where the operator Jg+n is defined by formula (TT4|) and the following convention is assumed: 
kg+i = 0. In case of initial states ([6]) that involve correlations such type (ITTl) cluster expansions 
(the kinetic cluster expansions) permits to take into consideration initial correlations in the 
kinetic equation. 

We give a few examples of recurrence relations flT7|) in terms of scattering cumulants flT3|) . 
Acting on both sides of equality flT7|) by the evolution operators ni=r^i('^' 0? ^^ obtain 

2ii(t,{y}) = QJi(t,{y}), 

s 

2l2(t, {Y},s + 1) = ^2{t, {Y},s + 1)+ QJi(t, {y}) J2 Mt, ii,s + l), 

il = l 

%sit,{Y},s + l,s + 2) = 5J3(t, {¥}, s + l,s + 2) + 
+2!5J2(t, {¥}, s + l)J2Mt,H,s + 2) + 



n=l 



+QJi(t,{F})(5^2t3(t,zi,s+l,s + 2) + 2! ^ ^2{t,zi,s + l%{t,i2,s + 2)), 

ii = l l=ii<i2 

where 2li+„(t) is the (1 + n)t^-order (n = 0, 1, 2) scattering cumulant (IT3|) . It is easy to prove 
that solutions of these relations are given by the following expansions 

s 

QJi(t, {¥}) = 2li(t, {¥}) = S,{-t, 1, . . . , s)3s{¥) J] S^{t, i), (18) 

i=l 
s 

^2{t, {¥}, s + l)= ^2{t, {¥},s + l)- 2ii(t, {¥}) Y, Mt, ^l,s + l), 
QJ3(t, {¥}, s + l,s + 2)= 2t3(t, {¥}, s + l,s + 2)- 2!2i2(t, {F}, s + 1) x 



^ 2t2(t, 21, s + 2) - 2li(t, {y}) ( ^ 2i3(t, ^1, s + 1, s + 2) - 

Ji=l il=l 

s s+1 s 

2!^ J]2t2(t,zi,s + l)2i2(t,i2,s + 2) + 2! ^ 2t2(t,zi, s + 1)212(^,^2, s + 2)). 



n=l«2=l l=n<42 

In the general case solutions of recurrence relations f lT7|) . i.e. the (1 + n)t/i-order (n > 0) 
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evolution operators Q3i+„(t), are defined by the expansions (IT2 



n— mi — ...— mj;_i 



QJi+„(t, {Y},X\Y)= n\ J2{-1)' E • • • E 7 ^ >< (19) 

x2ti+„_„i_..._„^(t, {y}, s + I, . . . , s + n - mi - . . . - ruk) X 

k ruj "'n—mi-...-m^+B—i s+n—mi~...—m,j 

11 Z^ • ■ • Z^ 11 TtJ _ U \j 

7 = 1 U—n U _n «i=l V"'n-mi-...-mj+s+l-ij "'n-mi-...-mj+s+2-i,- /• 

x2l,,.j _,i (f, i.-jS + n - mi - . . . - m,- + 1 + 

"'"'''s+n-mi-...-mj+2-ij' • • • j S + n — mi — . . . — m^ + fc^^.^.^^.^^ _„^^i_j^. J , 

s > 2, n > 0. 

This fact is verified as a result of the substitution of expressions (fT9|) into recurrence relations 
(TT7|) . It should be emphasized that the correlations generated by the dynamics of hard spheres 
and initial correlations connected with the forbidden configurations are completely governed by 
evolution operators ( IT9l) . 

Let us indicate some properties of evolution operators ( IT9l) . In the sense of the weak conver- 
gence of the space L^(R^'^ x ]R^^\Ws) the infinitesimal generator of first-order evolution operator 
flTSj) is defined by operator (jl]) 

limi(QJi(t,{F}) -/)/.= E ^ltihj)fs, 

where the operator C*^^{i,j) is defined by formula (j4]). 

In the general case, where n > 1, in the same sense as for evolution operators flTSl) it holds 
that 

iiniiQJi+„(t,{y},y\^)/,+„ = o. 

t^-o t 

Since in case of a system of non-interacting particles for scattering cumulants ( TT3|l the next 
equalities hold: ^i+nit,{Y},X\Y) = IlLi -^il"^' ^)^^(^) n?=i -^il^^iXo, n > 0, then we 
have 

s s 

^i+n{t,{Y},X\Y) = l[Si{-t,z)Js{Y)l[Si{t,j)5n,o, 

i=l j=l 

where Snfi is a Kronecker symbol. 

Similarly, at initial time t = it is true that: QJi+„(0, {Y} , X\Y) = 3s(Y)Snfl- 



Enskog Equation 11 

3.2 The derivation of the generahzed Enskog equation 

We establish that the one-particle distribution function defined by expansion (I8]),(|9]) in case of 
s = 1, i.e. 

oo ^ „ n+1 

Fi{t,xi) = y^ — / dx2...dxn+i^i+n{-t,l,---,n + l)Y\ F^{xi)X^3(i+u)\^^ , (20) 



n=0 



is governed by kinetic equation ( ITOj) or (fT5|) . 

In view of the validity of the following equalities for cumulants (^ of groups ([7]) for / G Lq 
in the sense of the norm convergence pQ , [1] 

liml2li(-t,l)/i(xi) = £t(l)/i(a^i), 

lim- / rfX22l2(-t, l,2)/2(Xi,X2) = / 6?X2£*nt(l, 2)/2(Xi, Xa), 

*^0 f^ JM3xR3 JR3xR3 

lim - / (ia;2 . . . dx^+i 2ti+„(-t, 1, . . . , n + l)/„+i = 0, n>2, 



*^0 ^ JlR3nx(R3n\W^) 



where the operators £^(1) and £*nt(l,2) are defined by formulas ([3]) and (jlj) respectively, as a 
result of the differentiation over the time variable of expression ( l20l) in the sense of pointwise 
convergence on the space L^(M^ x M^) we obtain [1], [H] 

-Fi(t,xi) = -(pi,— )Fi(t,xi)+ / rfx2Ct(l,2)5^-x (21) 

ot oil Jk3xm3 ^ n! 

„ n+2 

X / dX3...dx„+2 2ll+„(-t, {1,2},3, ...,n + 2)]TF{'(Xi)A:'IR;3(n+2)^w„+2■ 



i=l 



To represent the second term of the right-hand side in this equality in terms of one-particle 
distribution function fl20|) we expand cumulants ([9]) of groups ([7]) into kinetic cluster expansions 
f lT7|) for the case s = 2. Then we transform the series over the summation index n and the sum 
over the index ki to the two-fold series. As a result it holds 



n+2 



V - / dx^... dxn+2 2li+„(-t, {1, 2}, 3, . . . , n + 2) J„+2(l, . . . , n + 2) TT F^{xi) (22) 

^n!J(K3xR3)n f^l 

= V" — V" / rfxs . . . dx„+2+fci^i+n(i, {1, 2}, 3, . . . , n + 2) V . . . 

fcn+i ^ n+2 

• • • ,^5^„ ^n+2!(^n+l-fcn+2)!...(fcl-^2)! ^ 2il+'=".3-.-.„..-. (-^, ^, ^ + 3 + ^n+4-. 

n+2+fci 
. . . , n + 2 + fcn+3-j)'^l+fc„+3-i-fc„+4-i('?J) Q'n+3+fc„+4-i ) ■ ■ ■ ) Q'n+2+fc„+3_J J_ J_ -P^l 



^°(^.; 
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The last series can be expressed in terms of one-particle distribution function ( l20l) according 
to the following product formula for series ( l20l) 



n+2 oo fci fcn+i ^ 



„ n+2 

/ c^a;„_|_3 . . . dxn+2+ki I I 2li+fc„+3_i-fc„+4_i(— t, ^, n + 3 + A;„+4- 



n+2 

X 

'(M^xKS) . ^^]^ 

n+2+fci 

, n + 2 + /c„+3_j)A'i+fc„^3„-_fc„+4„i(gi, g„+3+fc„+4_i, . . . , gn+2+fc„+3_J J_ J_ -^] 



1 y^j, 



According to equalities fl22|) and fl23|) . and taking into account definition (j4]) of the operator 
£*jt(l,2) for t > 0, from equality ( 12T|) we finally derive 

— Fi(t,xi) = -(pi, — )Fi(t,a;i) + 

+0-^ V" — / dp2C??7 / (ixs . . . dxn+2{v, (Pi - P2)) x 



x{^^+4t,{r,2*^},3,...,n + 2)F,{t,qi,pl)F,{t,qi-ar],pl)l[Fi{t,: 

1=3 

n+2 

-5Ji+„(t, {1, 2+}, 3, . . . , n + 2)Fi(t, xi)Fi(t, gi + ar/,p2) JJ ^i(^' ^i)) ' 



«=3 



where we used notations accepted in equation (ITU|) . We treat the constructed identity for the 
one-particle distribution function as the kinetic equation for a hard sphere system. We refer to 
this evolution equation as the generalized Enskog kinetic equation. The collision integral series 
converges (see next section in the general case) under the condition: ||-Fi(t)||ii(]R3x]R3) < e~^. 

We emphasize that coefficients of the derived kinetic equation (ITU]) are determined by initial 
correlations connected with the forbidden configurations of particles. 

For quantum many-particle system the generalized kinetic equation has been formulated 
in [20] (see also a review pT]) and for systems of classical particles interacting via smooth 
potential it was stated by other method in [1], [22l|23] and for discrete velocity models in [21] . 

4 Marginal functionals of the state 

4.1 A mean value functional within the framework of the kinetic 
evolution 

By means of kinetic cluster expansions (ITTjl in case of s > 2, solution expansions ([8]) of the 
BBGKY hierarchy (Q for hard spheres can be represented in the form of the expansions with 
respect to one-particle distribution function O20p which is governed by the generalized Enskog 
equation (TTUl) . 
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Indeed, we expand cumulants (|9]) of groups ((71) from solution expansions (E]) into kinetic 
cluster expansions ( TT71) in case of s > 2. Then, transforming the series over the summation 
index n and the sum over the index ki to the two-fold series, we obtain 

s+n 



F,(t,xi,...,x,) = ^- / dxs+i...dxs+n'^i+n{-t,{Y},X\Y)3,+niX)Y[F^ 

oo ^ « oo ki 

= T- dxs+i . . . dxs+n 5Ji+„(t, {Y}, X \ y) V V . . . 

'^7l-\-S — l w « 

^'^Q ^n+sK^n+s-1 — kn+sV- ■ ■ ■ {h — A;2)! 7(R3xM3)fei 

n+s 

. . . dXn+s+ki J_ J_ ^l+kn+s+i-i-k„+s+2-i\'^^ Z, n + S + 1 + A;„+s+2-Ji . . . , 77, + S + Kn+s+l-i) 
i=l 
n+s+ki 

^ 111 {Xj)A,i^k„+s+i-i-k„+s+2-i\1i^ ?n+s+l+fc„+s+2-i) • • • 5 Q'n+s+fc„+s+i_J) 

i=i 

where we use notations accepted above. According to the validity of product formula fl23|) in 
case of the arbitrary index s > 2, in obtained expansion the series over the index ki can be 
expressed in terms of one-particle distribution function fl20|) . Hence the following equality holds 



oo ^ „ s+n 

y,-, dxs+i . . . dx,+„2li+„(-t, {Y}, X \ y) J,+„(X) n F{'( 

^ n\ J(M3xM3). fj^ 

oo -. „ s+n 

= y^- dXs+i...dXs+n'^i+n{t,{Y},X\Y)Y[Fi{t,Xi), 



where the {n + l)t/i-order evolution operator 5Ji+„(f) is defined by formula (IT2|) . 

Thus, if initial data is completely defined by a one-particle distribution function on allowed 
configurations, then the evolution of states governed by the BBGKY hierarchy ([1]) can be 
completely described by the generalized kinetic equation ( TTUj) and the sequence of marginal 
functionals of the state Fs(t,Xi, . . . ,Xs \ Fi(t)), s > 2, defined by expansions (TTB1) . In other 
words, for mentioned initial states all possible states of a hard sphere system at arbitrary 
moment of time can be described within the framework of a one-particle distribution function 
without any approximations. Therefore the statement of section 2.2 about an equivalence of 
stated Cauchy problems is valid. 

Typical properties for the kinetic description of the evolution of constructed marginal func- 
tionals of the state (TT6l) are induced by the properties of evolution operators (TT2|) . 

We indicate that expansions ( !T6|) of marginal functionals of the state are nonequilibrium 
analog of the Mayer-Ursell expansions over powers of the density of equilibrium distribution 
functions. In particular, if the one-particle distribution function is the Maxwellian distribu- 
tion function, then marginal functionals of the state fITB]) is transformed into the Mayer-Ursell 
expansions [T7], ^5\ f. 
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Owing that a sequence of solutions of the BBGKY hierarchy is represented in terms of a 
solution of the generalized Enskog equation by the marginal functionals of the state ( IT6l) . the 
average value of the s-ary marginal observable B^'^^ = (0, . . . , 0, as{xi, . . . , Xg), 0, . . .) [21] is 
defined in terms of the functional Fs{t,xi, . . . , Xs\Fi{t)), by the formula 

(S('*))(t) = — / dxi...dXsas{xi,...,Xs)Fs{t,xi,...,Xs\Fi{t)), 

and, in particular the average value of the additive marginal observable B^^^ = (0, ai{xi), 0, . . .) 
is determined by 

(S«)(t)= / dxMxi)F,{t,xi), (24) 



where the function Fi{t, xi) is a solution of the Cauchy problem (TT0|) -( 1TT|) . 

In fact functionals ( !T6|) characterize the correlations of states of a hard sphere system. We 
illustrate this statement by giving an example of the marginal functional from the collision 
integral of the generalized Enskog equation. In this case the marginal correlation functional 
G2[t,Xi,X2 I Fi{t)) is defined by the cluster expansion 

F2{t,Xi,X2 I Fi(f)) = Fi{t,Xi)Fi{t,X2) + G2{t,XuX2 I Fi(t)), 

and it is represented by the series similar to f ITB]) 

G2 {t, xi, X2 I Fi(t)) = (2li(t, {1, 2}) - I)F,{t, Xi)F,{t, X2) + (25) 

00 ^ „ n+2 

+ y.-i dx3 . . .dxn+2'^i+nit,{l,2},3, . . . ,n + 2)Y[ Fi{t,x,). 

^ nl 7{M3xR3)- fj[ 

Basically this cluster expansion gives the classification of all possible currently in use scaling 
limits [3], [6l[7] since it characterizes the deviation of the state from uncorrected state, i.e. the 
state satisfying a chaos property. In the scaling limits it is assumed that chaos property of 
initial state preserves in time, and hence the scaling limit means such limit of dimensionless 
parameters of a system in which the marginal correlation functional 6*2 (^) 1, 2 | Fi(t)) vanishes. 
According to definition f l25|) . it is possible, if particles of every finite particle cluster move 
without collisions. In consequence of the definition of a dispersion, for example, of an additive- 
type observable 

((5«-(i?«)(t))^)(t) = 

dxi{al{xi) - {B^^'^)'^{t))Fi{t,xi) + / dxidx2ai{xi)ai{x2)G2{t,xi,X2 \ Fi{t)), 



where the functional {B^^^){t) is determined by expression ([21]), the value of the dispersion 
is minimal for states characterized by marginal correlation functionals (1251) equal to zero, i.e. 
from macroscopic point of view the scaling evolution of hard sphere states is the evolution with 
minimal dispersion. 

Thus, in case of the Boltzmann-Grad scaling limit [T]-[6] the limit of solution (1201) of the 
generalized Enskog equation is governed by the Boltzmann equation for hard spheres and the 
limit of marginal functionals ( IT6l) are products of the limit one-particle distribution functions 
that means the propagation of a chaos. 
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4.2 On the convergence of a series for marginal functionals 

We prove the existence of marginal functionals of the state (fT6|l for Fi{t) G L^(R^ x M^). In 
case of smooth interaction potentials this fact was estabhshed in [22l[23] by other method. 

Owing to the fact that operator ([7]) is the isometric operator and the validity for |X\y| > 1 
the identity: ^p.(x\y)=M-Xi (""-'-) '^'~H|P| ~1)' = 0; ^^ cumulants ([9]) of groups ([7]) the following 
identity holds 

» s+n 

/ dXs+i...dXs+n^n{t,X\Y)Y[Fi{t,Xi) = 0, n>l. 



i=l 



As a result of the validity of these identities we can represent marginal functionals of the state 
flT6|) by the following renormalized series 

oo ^ „ s+n 

F,{t, xu...,xs\ Fi(t)) = V ^ / dxs+i . . . dxs+n^i+n{t, {Y}, X\Y)T\ F,{t, Xi),{2Q) 

where the (1 + n)th-oTdeT renormalized evolution operator QJi+„(t) is defined by the expansion 



n n-mi-...-m,fc_i 



%...(t,{r}.X\Y) = n'.^(-lf^... J: ^„^_'_ X (27) 

fc=0 mi=l ms,=l ^ -^ "^ 

X^l+n-mr~...~m^it, {Y}, S + 1, . . . , S + n - mi - ... -nik) X 
k mj A;^_i s ^ 

X n E • • • E n Tu — u — ^^i+^f..-.-^^..-. (^' ^^' 



=1 K'^s+l-i, ^s+2-ij)- ^ ^ 



s + n-mi-...-mj + l + A;^+2-i, , • • • , s + n - mi - . . . - m^ + kl^i_i^). 

and as above in (TT^ the next conventions are assumed: /c-J = mj, A;^_,_]^ = 0. 

To prove the convergence of the series of marginal functional fll6p we use the estimate of the 
cumulants 2li+„(t) of groups of operators ([7]) [18] 

s+n 
L} ■ 

1=1 



„ s+n 

/ dxi . . .cix,+n|2li+„(-t, {Y],X\ r) n Fi(t,x,)| < n!e"+lFi 

Jm3xM3)s+n 



In consequence of this estimate, for (1 + n)t/i-order renormalized evolution operator ( 1271) the 
inequality holds 

/s+n 
dxi . . .rfx,+„|^i+„(t, {Y},X\ y) n Fi(t,x,)| < n!||Fi(t)||^r x (28) 

n n n-mi-...-mfc_i fc m^ fc^_j ^ 



fc=Omi=l mfc=l i=l fe;'=0 fc^=0«j=l 
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So far as it is redefined that kl = rrij, kl^-^ = 0, we have 

s 

Then the following equality is true 

E • • • E 1 = T^^iyS'^^ + ^^^"^^ + 2)...{mj + s- 1). (30) 



fe^=o fci=o 



Indeed, using the equality (in case of m = it is the sum of the arithmetical progression) 

n ^ 

^ k{k + l) ...{k + m) = -n(n + 1) . . . (n + m + 1), (31) 

fc=i ^ 

as a result of the sequential computation of sums we obtain 

fc^=0 fc^=0 fc^=0 fcf_i=o fc^=0 fcs_2=0 

= • • • = (s-i)^ '^^ + ^)*^"^^' + 2) • • • ("^i + s - 1). 
According to estimates (l29l) and (l30ll . we have 

A" TTIt "* _ 1 s A" 

n E • • • E n e''=+i-'."''+^-'.-+' = n e^^+^s (m, + l)(m, + 2)...(m, + s-l) ^ ^g^) 

< TT r^i+'2s i^j + S — '^r ^ -TT 2mj+3s-l 

-AA (s-1)! ~ J-J- 

i=i ^ ^ i=i 

Thus, in consequence of estimate (l32l) inequality (128|) takes the form 



rfxi...rfx,+„|^i+„(t,{y},X\F)TTFi(t,Xi)| < (33) 

n n ?i-mi-...-mfe_i 

<n!||Fi(t)||^t"e"+^X]^^^^'"^^ E ••• E e"^^+-+™^ 



fc=0 m,i=l mfc=l 

Let US estimate the sums "^^1=1 ■ ■ ■J2^.=i~"~^''~^ ^"^^~^"'~^'^''- Calculating the sum over 
index ruk as the sum of a geometric progression, we obtain 

n n-mi-...-m,fc_i n n-mi-...-mfc_2 

Y, ■■■ E e'"^+-+'"* < e"+i J2 ... Yl 1- 

■mi=l m;j=l mi = l 'mfc_i=l 
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Then, applying equality (13T1) . as a result of the sequential computation of sums the following 
inequality holds 

n n-mi-...~mk-i n n-mi--...-mk-z 

^ . . . ^ gmi + ...+mfe < e^+i ^ . . . ^ (^ _ m^ _ . . . _ ^^2) = 

mi=l mfe=l m,i=l mfe_2=l 

^ , (n - A; + 2) ■ ■ ■ (n - l)n 
(fc-1)! 

According to the last inequality, for the majorant of estimate (1331) we get 

n n n-mi-...-m.fe_i 

fc=0 r)ii=l mfc=l 

< 1 + e-+3s + gn+i y- „fc(3s-i) (^-fe + 2)---(n-l)n ^ 
~ ^^ (A;-l)! ~ 



fe=2 



< 1 + e"+3- + e"+i V e'^^^^^-i)— < 1 + e^^+i V 

~ ^^-^ (A;- 1)! ~ ^^-^ 

fc=2 ^ ^ A;=l 



gfe(3s-l)_ 



Thus, in consequence of inequalities (15^ and (151|1 for marginal functionals of the state 
the following estimate holds 

00 00 n 

\\Fs{t I Fi(t))||^, < Y, ||Fi(t)||;|V+2 + 5] ||Fi(t)||;t\«^e'=(=^^-i) < (35) 

n=0 n=l fc=l 

00 00 

^ \\W (+\\\'' p2\^ II p /.Nll'^ n I II p /.N||« p3s+2 Y^ II p /.Njj'i ^(3s+2)n 

^ IriWlLie 2^Piltj||iie +||i^i(t)||^ie 2^|KiW|Lie^ ^• 



Hence, the sequence of marginal functionals of the state Fgit^xi, . . . ,Xs \ Fi(t)), s > 2, exists 
and represents by converged series ( TT6l) provided that 

||i^i(t)|Ui(R3xR3)<e-(3^+2)_ (3g) 

On the other hand, owing the estimate for series ( l20l) 

00 
||^i(^)||li(ir3xir3) < ll^i°IUi(iK^xiR3)e^^e"'||Ff||2i(iR3xK3), (37) 

n=0 

and according to estimates fl36|) and fl37j) . we conclude that the Enskog collision integral exists 
under the following condition on initial data 

g-10 

11^1 ||l1(M3xK3) < _Q . (38) 

-L ~r t- 
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5 Some properties of the collision integral 

5.1 The Markovian generalized Enskog equation 

To state the Markovian generalized Enskog equation we first represent the coUision integral 
1-GEE of the generalized Enskog equation (TTO!) as an expansion with respect to the Boltzman- 
Enskog collision integral (in case of t > 0) 

'^BEE = 0- / dp2dr]{r],{pi - P2)){Fiit,qi,pl)Fi{t,qi - ar],pl) - Fi{t,Xi)Fi{qi + ar],p2)), 

where the momenta pl,P2 are defined by equalities 1^. We observe that such expansion of the 
collision integral Xqee is given in terms of marginal correlation functionals fl25|l . Indeed, in 
view of that it holds: (2ii(t, {1«,2^}) - I)Fi{t,qi,p\)Fi{t,qi±ari,p\) = 0, we have 



oo 



'^GEE — ^BEE + 7 ,'^GEE = (39) 

n=l 

= Ibee + (^^y2~i / ^^2C^^ / dx3... dXn+2 iv, {Pi -P2)) X 

n+2 

x{^,+4t,{l*,2*_},3, . . . ,n + 2)F,{t,q,,pl)Fi{t,q, - ar],p*^)l[F,{t,Xi) - 

i=3 

n+2 

-^i+n{t,{l,2+},3,. . .,n + 2)F^{t,x,)F^{t,q, + ar],p2)Y[F,{t,Xi)), 



«=3 



where notations accepted in equation (fTOj) are used. 

On the kinetic (macroscopic) scale the typical length for the kinetic phenomena is the mean 
free path. Then, observing that in the kinetic scale of the variation of variables [2], [7] the 
groups of operators ([7]) of finitely many hard spheres describe a fast evolutionary processes and 
hence they depend on microscopic time variable e~^t, where e > is a scale parameter (the 
ratio of the collision time to the mean free time), the dimensionless marginal functionals of the 
state are represented in the form: Fs(^e~^t, xi, . . . ,Xs \ Fi{t)), s >2. In the formal limit e — )■ 0, 
(the Markovian approximation) the limit marginal functional of the state Fs{xi^ ■ ■ ■ ,Xs \ Fi{t)) 
is represented by expansion ( 1T6l) (and also ( l25l) ) with the limit generating evolution operators 
\im.e^o^i+n{£~^i) , n>0, for example, 

lim QJi {e-H, {Y}) = lim §1 (e^^t, {Y}). 

We note that the Markovian limit of the first two terms of expansions flT6|) of dimensionless 
marginal functionals of the state coincide with corresponding terms of particular solution of the 
BBGKY hierarchy constructed by N.N. Bogolyubov ^, [T7], [22] on basis of the perturbation 
method with the use of the condition of weakening of correlation. 
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In the formal Markovian limit e — )■ from the generalized Enskog equation ( ITOl) we derive 
the Markovian generalized Enskog equation 

— F^{t,Xi) = -{pi,—)F^{t,X,)+lBEE+ (40) 

+0-M dp2dr]{r],{pi-p2))y2—i (ixg . . . rfa;„+2(2Ji+n({l*, 21}, 3, . . . 

n+2 

. . . ,n + 2)F,{t,qi,pl)Fi{t,q, - cxr],p;)l[F,{t,x,) - 

1=3 

n+2 

-QJi+„({l, 2+}, 3, . . . , n + 2)Fi(t, Xi)Fi(t, q^ + ar],p2) JJ Fi{t, Xi)) , 



i=3 



where by QJi+ndl**, 2!|_}, 3, . . . , n + 2) we denote as above the limit generating evolution opera- 
tors: QJi+„ = \mv^^Q^i+n{e~H), n > 0. 

Thus, in the kinetic scale of the variation of variables the structure of the collision integral 
of the generalized Enskog equation (HOj) describes all possible correlations which are created 
by hard sphere dynamics and by the propagation of initial correlations stipulated by forbidden 
configurations. 

We remark that for the first time two terms of the Markovian collision integral in case of 
smooth interaction potentials was established by N.N. Bogolyubov [I?] within the framework 
of the perturbation theory and using an analogy with the virial equilibrium expansions, these 
terms of the collision integral in non-perturbative form were formulated by M.S. Green and 
R.A. Piccirelh [26], [28] and by E.G.D. Cohen [29]. 

To derive the Bogolyubov collision integral for hard spheres from (140!) we apply analogs of 
the Duhamel equations [27j to generating operators of marginal correlation functionals fl25p . 
Indeed, if /s G L^{Mj^ x {W \ W3)), for the scattering cumulant 212(5-^, 1, 2, 3) ([I3D an analog 
of the Duhamel equation holds [1| 

^2{e-H, {1,2}, 3)/3(xi, xs, X3) = / drS2{-T, 1, 2)5i(-r, 3) V £^(^1, 3) x 

-^0 n=l 

3 

xS'3(-£"4 + r, 1, 2, 3)A']R9\W3 W Si{e'H, i2)h{xi., X2, X3), 

i2=l 

and, consequently, for the first term of the collision integral expansion we obtain 

-^GEE ~ -^GEE,+ -^GEE-^ 
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where we denote by Iq^e ± ^^^ following parts of the collision integral Iqee 

Igee,± = <^' [ dp2dv{v,{pi-P2)) [ dxs r dTS2i-r,l\2i){{Cl,{lKs) + 
JR^xsl Jr^xr^ Jo 

+Ct(2L 3))S3{-e-H + r, 1«, 2^, 3)X2{qu 93)^2(91 ± ^r/, ga) x 

xSiie'H - T, l^)Si{e-H - T, 2i)Si{e-H - r, 3) - S2i-e-H + r, 1«, 2i)Si{e-H - r, 1«) x 

xS,ie-H - r, 2i)(£*,,(l«, 3)52(-e-'t + r, 1«, 3)A'2(gi, g3)5i(e-H - r, l»)5i(e-H - r, 3) 

+Ct(2L 3)52(-e-'t + r, 2^ 3)A'2(gi ± ar/, g3)5i(e~H - r, 2^) x 

xSiie-H-r,3)))S^ir,l^)Si{r,2i)S,iT,3)F,it,qup\)F,it,qi±av,pl)F,it,Xs). 

In the Markovian limit e — )■ this expression coincides with the corresponding correction to 
the Boltzmann-Enskog collision integral constructed by the perturbation method with the use 
of the weakening of correlation condition in [17], [25]. 

5.2 To the classification of the Enskog- type kinetic equations 

We establish links of the Markovian generalized Enskog equation ( l4Up with the revised Enskog 
equation [SHTU]. The classification and links of the Enskog-type kinetic equations derived upon 
some heuristic arguments is given in [8], [T2] . 

The collision integral of the revised Enskog equation has the form [8HTO] 

^REE = crM dp2d'r]{r], {pi - P2)){g2{qiAi - (^V I Fi{t))Fi{t,qi,p\)Fi{t,qi - ar],pl) - 

jR^xSl 

92iqi,qi + crr] I Fi{t))Fi{t,qi,pi)Fi{t,qi +r]a,p2)), 
where the functional g2{qi, ^2 | Fi{t)) is defined in terms of the formal Mayer cluster expansion 
g2iquq2\F,it))= (41) 

00 ^ „ fc+2 

= X2iqi, ^2) (1 + y] 77 / dxs... dxk+2Vi+ki{qi, 92}, ^3, • • • , qk+2) TT Fi{t, Xi)) , 

and Vi+ki{qi, ^2}, qs, ■ ■ ■ , qk+2) is the sum of all graphs of k labeled points which are biconnected 
when the Mayer function f{qi, ^2) = ^^2(^1, g2) — 1 is added, ^'2(^1, q'2) = '^rb\W2 is the Heaviside 
step function of allowed configurations of two hard spheres. Taking into account that the 
Boltzmann-Enskog collision integral is defined by the first term of expansion (SI]), we denote 
the collision integral of the revised Enskog equation in the form of the corresponding expansion 

00 

^REE = ^BEE + / S^REE,+ ~ ^REE-)^ (42) 

fe=l 

where 

^REE,± = ^ / dp2dr]{r], {pi - P2)) I dqs . . . dqk+2V{{qi, qi ± ar]},q3, . . . , qk+2) x 

jR^xSl JR^k 

k+2 

xFi{t, qi,p\)Fi{t, qi ± ar],pl) J] Fi(t, ; 



k+2 
i=2 
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We compare term by term collision integral (H2l) of the revised Enskog equation and the 
collision integral ( 139|) of the Markovian generalized Enskog equation ( 140|) by the example of the 
first terms of these expansions. In case of collision integral ( H2|) this term has the form 

^REE,± = (^^ dP2dv{v,{Pi-P2)) dx3fiqi,q3)f{qi±crri,q3) X (43) 

xFi{t,qi,p\)Fi{t,qi±ar],pl)Fi{t,X3) = 

= crM dp2dr]{r], (pi - ^2)) / dx3{X2{qi, q3)X2{qi ± ar], gs) - 

-^2iquq3) - X2{qi ±0-77,^3) + l)Fi{t,qi,p\)Fi{t,qi ± crr],pl)Fi{t,X3), 
and in case of collision integral fl39|) the corresponding term has the following form 

^gIe± = (^' [ dp2dr]{7],{pi-p2)) [ dx3\im(^2{£-H,{lK2i},3) - 

-2ti(e-4, {1«, 2i})2t2(£-4, 1«, 3) - 2ti(£-4, {!«, 2^})2t2(£-^t, 2^, 3)) x (44) 

xFi{t,q,,p{)F^{t,q,±a7],pl)F,{t,X3), 

where we use notations accepted in equation (fTOj) and the evolution operator 2l„(e~^t) is a 
corresponding order scattering cumulant (TT3|) . 

We transform expression (1H|) to the form adopted to the structure of expression 043p . With 
that end in view we introduce the notion of the scattering operator 

n 

Snil, . . .,n) = limS'„(-£:"^t, 1,. . . ,n)Y\ Si{6~'^t,i), n > 1, 

where the evolution operator Sn{—£~^t) is defined by ([7]). Then, according to definition (fT3!) of 
scattering cumulants and definition (jTj), for the first summand in expression (jH]) the following 
equality holds 



Iim2l2(£-'t, {!«, 2^}, 3) = [fiqi, q3)f{qi ± ^77, ^3) + /(^i, ^3) + f{qi ± cjt], ^3) + 
+1)^3(1«, 2L 3) - lim S2{-e~\ l\ 2i)5i(-£-4, 3) (/(gi, g3)/(?i ± ar^, ^3) + 
+/te, 93) + f{qi ± ^^, gs) + 1)S^{e-H, 1»)5i(£-4, 2^)5i(£-4, 3), 



where we take into consideration that S3{e ^t, 1^, 2^., 3) {f{qi, q3)f{qi±ari, q3) + f{qi, q3) + f{qi± 
'^Vi Q3)) = 0- For other summands in expression (jH]) it holds, respectively 



]im^^{E-H,{lK2i})^2{e~H,lKS) = S2{lK2i){f{q,,q3) + l)S2{lKS) - 

lim S2i-e-H, 1«, 2^ ^i(-£-it, 3) (/(gi, ^3) + l)5i(£-4, l^)Sii6~h, 2i)Si{e~H, 3), 



£^•0 

-h 

£-S>0 



lim2ti(£-it, {!«, 2^})2l2(£-it, 2i, 3) = 52(1«, 2»±) (/(gi ± ar^, 53) + 1)52(2^, 3) - 

- lim S2{-e-h, 1«, 2i)Si{-e~H, 3) (/(gi ± ar/, ^3) + l)Si{e~H, l^)Si{e'H, 2i)Si{e~H, 3). 
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As a result of the validity of these equalities, taking into account that 5*2(1'', 2?|_) = / is the 
identity operator and, owing to the following equality 

/ dxs lim S2i-s-\ l\ 2i)S,i-e-\ 3)/(gi, q^)fiqi ± ar/, ^3) x 

xS,{e~H,li)Si{e-H,2i)S,{e-H,3)F,{t,q^,p{)F,{t,q^±ar],pl)F,{t,X3) = 0, 
expression (jH]) takes the form 

^GEE,± = ^ / dp2dr]{r], (pi - pa)) / dx3{X2{qi,q3)X2{qi ± (tt], q^) x (45) 

xS3{lK2i,3)-X2{qi,q3)S2{l\3)-X2{qi±cxv,q3)S2{2i,3) + l)x 
xFi{t,qi,p\)Fi{t,qi±ar],pl)Fi{t,X3). 

It is clear that last expression 0451) coincides with corresponding term (1431) of the collision 
integral expansion of the revised Enskog equation, if we neglect the collisions which occur in 
the process of the evolution of hard spheres. 

As we note above in section 4.1 expansions of marginal functionals of the state flT6|) are 
nonequilibrium analog of the Mayer- Ursell expansions over powers of the density of equilibrium 
distribution functions [TTj, [2^126]. Hence, the structures of expansions of collision integral (1421) 
of the revised Enskog equation and collision integral (l39l) of the Markovian generalized Enskog 
equation (HO!) are the same, so observations given for the first terms of these collision integral 
expansions take place for every term. In fact every term of collision integral (139!) give us the 
classification of correlations generated by pair collisions of hard spheres. 

Thus, in case of neglecting of collisions in particle clusters of more than two hard spheres, 
the collision integral of the generalized Enskog kinetic equation fHOl) transforms to collision 
integral (H2l) of the revised Enskog equation. 

6 The Cauchy problem of the generalized Enskog equa- 
tion 

6.1 An existence theorem 

We consider the abstract Cauchy problem flTU]) - flTT]) ( flT^ . flTTl) ) in the space of integrable func- 
tions L^(M'^ X M.^). The following statement is true. 

Theorem 1. A global in time solution of the Cauchy problem of the generalized Enskog equation 
l[T^) - l[Il\) ( in^) . ni]) ) is determined by the expansion 

00 _ „ n+l 

Fi{t, xi) = V" — / dx2... dxn+i 2li+„(-t, 1, . . . , n -h 1) TT F^{xi)Xj^3(i+n)\^^ (46) 

where the cumulants 2li+„(— t), n > 0, of groups of operators ^ are defined by formula ^. If 
||-^?||li(r3xR3) < e^^^{l + e^^)^"'^, then for F° G -/^o(M'^ x M.^) it is a strong solution and for an 
arbitrary initial data F^ G L-'^(]R'^ x M'^) it is a weak solution. 
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Proof. Let F° G Ll{R^ x R^). Series ( H6l) converges in the norm of the space L^(R^ x M^) 
provided that ( 138|) . We use the result of the differentiation over the time variable of expression 
(146|) in the sense of pointwise convergence on the space L^(R^ x R^) from section 3.3. Then for 
F^ e LJ(R3 X R3) it holds that 

lim / dxi|i-(Fi(t + At,xi)-Fi(t,xi))- (47) 

oqi ^ n! yiR3xR3 

^ n+2 

X / C/Xa . . . dXn+2'^l+n{t, {1, 2}, 3, . . . , 72 + 2) TT Fi(t, X,)) I = 0, 

7(R3xR3)" 7Li 



where the operator C*^^{1,2) is defined by formula @ and abridged notations are used: the 
symbols Fi(t, xi) and 11^=1 -^i(^; ^i) ^^^ implied series fH6|) and for s = 2 series (!23|) . respectively. 

The validity of equality f l47p means that a strong solution of initial- value problem fl 10 1) - fill I) 
is given by series ( H^ . 

Let us prove that in case of initial data F^ G L^(R'^ x R'^) a weak solution of the initial- 
value problem of the generalized Enskog equation fITOl) is represented by series fH6l) . With this 
purpose we introduce the functional 

((^i,Fi(t))= /" dxm{xl)F^{t,x^), (48) 

Jr^xM.3 

where ipi G i^o(IR'^ x I'^'^) is continuously differentiable bounded function with compact support 
and the function Fi{t) is defined by series (1461) which is convergent series in the norm of the 
space L^(R^ x R^) under condition (1381) . This functional exists since the function Fi{t) is 
integrable and (pi is a bounded function. 

Using expansion P6|) . we transform functional fHHj) as follows 



{^i,F,{t)) = 

n+l 

(l + n) 



V'-T / dxi...rfz„+iv?i(xi)2ti4.„(-t, l,...,?2 + l)TTF°(2;i)A'i 

^ '^l J(R3xIR3)n+l fj^ 

oo ^ „ n+l 

V] — / dxi . . . dx„+i2ti+„(t, 1, . . . , n + l)(/5i(xi) TT F°(xi)A:']K3{i 

^'^l J(R3xR3)n+i f^ 



where the cumulant 2li+„(t), n > is adjoint to the cumulant 2li+„(— t), n > 0, in the sense 
of functional dH]). For F° G L^(R3 x R^) and ipi G Lo(R^ x R^) in the sense of the *-weak 
convergence on the space L(R^ x R'^) the following equality holds 

oo ^ 1 

lim 52 ^ / dxi... dxn+i — (2ti+„(t + At) - (49) 

n+l 

-2ti+„(t))(/)i(2;i) JJ i^i°(a;i);fR3(i+n)\Wi+„ = 

i=l 
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°° 1 /" 

= {Ciipi, Fi{t)) + V" -r / dXi... dXn+2'Cint{l, 2)(pi{xi) X 

^ "-! J(IR3xIR3)"+2 

n+2 

x2li+fe(-t, {1, 2}, 3, . . . , n + 2) J] F{'(x,)A'j,3(„+2)\w„+,- 



j=i 



In obtained expression (H9l) the operator £int(l, 2) is adjoint to the operator >Cj*jjt(li 2) and it is 
defined on Lo{R^ x (M^ \ Wa)) by the formula 



/ dxidx2Cmtil,2)(p2{xi,X2) = (t"^ / dxi / 

J(R3xIR3)2 jR3xR3 JlR3 



dp2dr]{r],{pi-p2)) x (50) 



For F}* G L^(]R^ x R^) and (fi G Lo(M^ x R^) the hmit functional (HHD exist. Using equahty 
fl23|l . we transform the second functional in (H9l) to the form 



1 /■ 
V - / dxi . . . (ia;„+2Ant(l, 2)vJi(xi)2li+fc(-t, {1, 2}, 3, . . . , n + 2) x (51) 

„^;^ "-' J(R3xR3)"+2 

n+2 oo ^ „ 

^ TT ^l'(^i)'^R3("+2)\W„+2 = y^ ~r / rfXi(iX2£int(l, 2)v?l(Xi) X 

r=l ^ '^^ i(R3xR3)2 

» n+2 

X / rfx3...(ia;„+25Ji+„(t, {1,2},3, ...,n + 2) TTFi(t,Xi). 

7(R3xR3)n ^^-^ 

Therefore, in consequence of equalities ( H^ .l lFT]) . for functional ( HSjl we obtain the equality 



n+2 

da;3...rfx„+22Ji+„(t) 

n=0 '*• -^(K^xR3)" .^^ 

or in the explicit form 



» oo _. „ n+2 

/ dxidx2Cmt{l,2)ipi{xi)y2~i (ia;3...rfx„+22Ji+„(t) TTFi(t,z) 

7m3xR3)2 ~:; ^! J(R3xR3)" ,J7 



1 r Pi 

-((/9i,Fi(t))= / rfxi(pi,^^)(/Pi(xi)Fi(t,xi)+ (52) 

at jR3xR3 OQi 

+a^ dxi dp2dr]{r],pi-p2){ipi{qi,pl) -(pi{qi,pi))y2~; dx^ . . . 

JR3x]R3 Jr3xS2_ ^^^ nl J(K3xR3)n 

n+2 

. . . dxn+2'^i+n{t, {1, 2+}, 3, . . . , n + 2)Fi(t, xi)Fi(t, gi + ar],p2) JJ Fi(t, x,). 



j=3 



This equality ( l52l) means that series (l46l) for arbitrary initial data F° G L^(]R'^ x M^) is a weak 
solution of the Cauchy problem of the generalized Enskog equation (iTOi) - (ITIi) . 

D 
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6.2 A weak solution in the extended sense 

For the Cauchy problem (TT0|) - (TTT1) we introduced the notion of a weak solution in the extended 
sense. Let ((/?, F(t | Fi(t))) be the functional 

(v3,F(t I Fi(t))) = V'-r / dxi . . .dxs(Ps{xi, . . . ,Xs)Fs{t,xi, . . . ,Xs \ Fi{t)), (53) 

where (p = {(po, ipi{xi), . . . , (ps{xi, . . . , Xs), ■ ■ ■) is a finite sequence of bounded infinitely times 
differentiable functions with compact supports (ps £ ivo(K^'' x (R^'^ \ W^)) G L{W'^ x (M^"* \ 
Ws)), s > 0, and elements of the sequence F{t, \ Fi{t)) = (l, Fi(t, xi), F2{t, xi, X2 \ Fi{t)), . . . , 
Fs{t,xi, . . . ,Xs I Fi(t) ),...) are defined by formulas P6|) and flT6l) for the first and other 
elements respectively. Owing to estimates fl37|) and fl35|) . this functional exists provided that: 

On functions ips ^ i^o(J^^'' >< (J^^* \ Ws)), s > 1, it is defined the operator B^ which is adjoint 
to the generator of the BBGKY hierarchy of hard spheres ([1]) for t > 0, 

s 
{B''ip)s{Xi,...,Xs) = CsiY)(p!s{Xl,...,Xs) + ^ CintUlJ2)Vs-l{ixi, . . . ,Xs)\iXj^)), 

where the operator £int(ji, J2) is defined by ( !50|) and the Liouville operator Cs(Y) is adjoint to 
operator ([3]) in the sense of functional ( 153|) [Ij . 

It is said to be that expansion f H6l) is a weak solution of the Cauchy problem (ITOi) -( TTT|) in 
the extended sense, if for functional (153!) the following equality is true 

j^{^,F{t\F,{t))) = {B^^,F{t\F,{t))). (54) 

In view of estimates (!35|) the functional on the right-hand side of this equation exists provided 
that (|36|). 

To verify this definition we transform functional (l53l) as follows [21] 



(^,F(t|Fi(t))) = 

00 ^ „ s ^ s 

^=0 >^(R xK ) n=0 ^ ^ ji^...5^j,_„=l Zcy\(ji,...,i.-„) 

s 

x5',_„+|z|(t, (ji, . . .,js-n) U Z)ip,_ri{Xj^, . . . ,Xj^_„) JjF°(Xi)A'iR3.\W,, 



i=l 



where X]zcr\(?i ?«_„) ^^ ^ ^^'^ °^^^ ^^^ subsets Z C ^\(ii, • • • , js-n) of the set Y\{ji, . . . , js-n) C 
Y = (1, . . . , s). For "f^ e L\M.^ x R^) and y?, G Lo(M^" x (M^" \ W,)), s > 1, this functional 
exists. 

Skipping the details, as a result for ips E Lo(R^'^ x (M^* \ Wg)), s > 1, the derivative of 
functional f l53|) over the time variable in the sense of the *-weak convergence on the space of 
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sequences of bounded functions takes the form 

d °° I f 

-{ip,F{t\Fi{t)))=y^- dxi...dxs{Cs{Y)^s{xi,...,Xs)+ (55) 

"■t ^ ■S! 7(R3xR3)s 

s 
+ 5Z ^mt{jl,J2)Vs-li{Xu...,Xs)\iXj^)))Fs{t,Xi,...,Xs \Fi{t)). 

According to the notion of a weak solution in the extended sense fl54p . equahty fl55p means 
that for arbitrary initial data F^ G L^(R^ x R^) a weak solution of initial- value problem of the 
generalized Enskog equation f lTOj) is determined by series fH6|) . In fact, we have also proved 
that marginal functionals of the state flT6l) are weak solutions of the BBGKY hierarchy of hard 
spheres. 

7 Conclusion 

Thus, in case of initial data (|6]) solution (|8]) of the Cauchy problem (II])-© of the BBGKY 
hierarchy for hard spheres and a solution of the Cauchy problem of the generalized Enskog 
equation flT(I]) -f lTT]) together with marginal functionals of the state fITB]) give two equivalent 
approaches to the description of the evolution of a hard sphere system. We emphasize that in 
addition, the coefficients of the generalized Enskog equation ( ITOj) and marginal functionals of 
the state (TT6|) are determined by the operators of initial correlations specified by the forbidden 
configurations of hard spheres. 

The structure of the collision integral expansion of the generalized Enskog equation (ITOl) is 
such that the first term of this expansion is the Boltzman- Enskog collision integral and the next 
terms describe all possible correlations which are created by hard sphere dynamics and by the 
propagation of initial correlations connected with the forbidden configurations. 

On the basis of the derived generalized Enskog equation flTOj) we state the Markovian gener- 
alized Enskog equation (HU]) . In case of absence of correlations which generated by hard sphere 
dynamics the collision integral of this kinetic equation transforms to the collision integral of the 
revised Enskog equation [8]. Therefore the kinetic evolution of hard spheres described by the 
revised Enskog equation makes it possible to take into consideration only the initial correlations 
specified by the forbidden configurations of particles. 

In the end it should be emphasized that the kinetic evolution is an inherent property of 
infinite-particle systems. In spite of the fact that in terms of a one-particle marginal distribution 
function from the space of integrable functions can be described a hard sphere system with the 
finite average number of particles, the Enskog kinetic equation has been derived on the basis 
of the formalism of nonequilibrium grand canonical ensemble since its framework is adopted to 
the description of infinite-particle systems in suitable functional spaces [T] as well. 

We remark also that developed approach is also related to the problem of a rigorous derivation 
of the non-Markovian kinetic-type equations from underlaying many-particle dynamics which 
make possible to describe the memory effects of particle and energy transport, for example, the 
anomalous transport in turbulent plasma, the Brownian motion of macroparticles in complex 
fiuids. 
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